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We investigate the microwave absorption arising from inhomogeneity in the superfluid density of a model high- 
Tc superconductor. Such inhomogeneities may arise from a wide variety of sources, including quenched random 
disorder and static charge density waves such as stripes. We show that both mechanisms will inevitably produce 
additional absorption at finite frequencies. We present simple model calculations for this extra absorption, and 
discuss applications to other transport properties in high-Tc materials. Finally, we discuss the connection of these 
predictions to recent measurements by Corson et. al.^ of absorption by the high-temperature superconductor 
Bi2Sr2CaCu2 08+i in the THz frequency regime. 



1. Introduction 

The high-Tc superconductor 

Bi2Sr2CaCu208+5 shows a remarkably strong ab- 
sorption in the microwave regime, even far below 
Tc0,^. In conventional, low-Tc, s-wave super- 
conductors, there is no such background, because 
there is no absorption below the energy gap for 
pair excitations, 2A. But in high-Tc materials, 
which are thought to have a (1^2 _y2 order param- 
eter, the gap vanishes in certain nodal k direc- 
tions. Hence, gapless nodal quasiparticles can be 
excited, and hence can absorb microwave radia- 
tion, at arbitrarily low temperatures. However, 
experiment suggests that this absorption is much 
stronger than expected from the quasiparticles 
alone HJ. 

In this paper, we show that such extra absorp- 
tion can be produced if the superfluid density 
n(x) within the ab plane is a function of posi- 
tion X. Such spatial variation can be produced, 
e. g., by quenched disorder, or by charge den- 
sity waves (e. g. stripes). The extra absorption 
is most likely in the microwave frequency range. 



A similar absorption was found by us in a pre- 
vious paper Q, in which the inhomogeneity was 
described as static fluctuations in the Josephson 
coupling between superconducting grains. 

2. Microwave Absorption in an Inhomoge- 
neous Two-Dimensional Superfluid 

We consider a two-fluid model of a supercon- 
ductor in two dimensions, with local conductivity 



(1) 



where q — 2|e| and m* is twice the electron 
mass TOg. n<,(x, oj) is assumed spatially varying, 
while Gqp is taken to be spatially uniform. We 
wish to calculate the complex effective conduc- 
tivity (Te(w). Gqp might be a contribution from 
the nodal quasiparticles, while Ug represents the 
perfect-conductivity response of the superconduc- 
tor. It can be spatially varying because of the 
very short in-plane coherence length. 

We employ the Kramers-Kronig relations sat- 
isfied by cre(ijj), namely 
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O-el(w) = /o°° ^ dJ + (Too, and CTe2 = 

_ r° Zsii^O^rf^', where P means 
"principal part of," Ue — cFei + *fe2, and Us^e is 



the effective superfluid density [the value of 
can be found from the residue of ae{uj) at w = 0]. 
At very large frequencies a;, the equation for CTe2 
becomes 0^2 ^ ^^-^ + ;|j S^[(Jei{uj') - a^]duj' . 

2.1. Frequency- Independent (jqp 

If o'gp(cLi) is real and frequency- independent, 
then (Too = o'qp. Then at high frequencies, the lo- 
cal complex conductivity, aoo + i<f'ns{yi)/(jn*uj), 
has only small spatial fluctuations. Thcng 

1 {{5af) 

6a{uj,x) EE CT(a;,x) - CTat,(w), where crQ.u(tj) = 
CToo + iq^ns^avl{m*uj), (5ct(cj,x) = cr(a;,x) - 
CTa^(w), and (...) is a space average. Since only 
Us, and not cr^p, is fluctuating, this expression 
simplifies to 



^ 1 {q^ / 771* f{{5nsf) 

m*uj 2 w^fjoo + iq'^ums^av/in* 



.(3) 



At large lu, the imaginary part of cTg to leading or- 
der in 1/uj is simply cre2 ~ q'^ris^av/ (rrfui). Equat- 
ing this expression to the right-hand side of the 
Kramers-Kronig expression for (7^2 , we finally ob- 
tain 

/ [a el {oj' )- (T 00] duj' = {Us^av - ris^e) (4) 

Jo 2m* 

If Tig (x) is spatially varying, the right-hand side is 
always positive, whence there will be an additional 
contribution to cjeii'-^), beyond Uqp. 

At small uj, eq. (||) implies ae2 ~ ' ' 
and thus n,,,. 



{{5nsf)/{2n 



^ ) 

s.av I 



{oeiiuj') ~ (J^\dJ 
4m 



Then eq. becomes 



(5) 



For fixed (((5ns)^)/"-s am this integral is propor- 
tional to ns,a«- That is, the extra integrated fluc- 
tuation contribution to Uei , is proportional to the 
average superfluid density. A similar result has 
been reported in experiments |]l|,^. 



2.2. Drude aqp{uj) 

For a Drude aqp{ijj) = cro/(l — iujr), we can 
carry out a similar analysis]^. In the case of weak 
fluctuations in n,, the result islpl 



[ae,i{io') ~ aqp^{uj')]duj' 
^^-{{qVm*f{{5nsf)) 



(6) 



q^ns^av/m* q^Us^av/m* + Uq/i 

In the limit ctq/''' 3> q^Us^av/fn* , this expres- 
sion reduces to the right-hand side of eq. (||). 
Thus, in this regime, the extra spectral weight 
is indeed proportional to the average superfluid 
density n^^av 

2.3. Tensor 

Next, we consider a tensor superfluid density, 
as would be expected in a superconducting layer 
containing quenched charge density waves such 
as a charge stripes 0. In this case, the (2 x 2) 
superfluid density tensor should have the form 
n"'^(x) = _R^^(x)nf _R(x), where is a diago- 
nal 2x2 matrix with diagonal components Us^Ai 
Us^B, and -R(x) is a position-dependent 2x2 rota- 
tion matrix, describing the relative orientation of 
the charge density wave or stripes. If the stripes 
have either of two orientations along the crystal 
axes of the layer, with equal probability, then the 
effective superfluid density n^^e will be a scalar. 

The arguments of the previous subsections can 
readily be transferred to the tensor case. For 
a two-fluid model with a frequency-independent 
quasiparticle conductivity which has the same 
value (Too in both the A and B directions, the 
effective scalar conductivity fTe(w) again satisfies 
Kramers-Kronig relations, and one again obtains 
the sum rule (]^, where Ug^av is now the rota- 
tional average of a diagonal element of ng-afi- If 
the principal axes of Ug-.a/i point with equal prob- 
ability along the two symmetry directions of the 
Cu02 plane, as mentioned above, Ug^av = ins,A + 
ns,B)/2. Likewise, if the principal axes were to 
point in any direction in the plane with equal 
probability (a circumstance which seems unlikely 
for a stripe phase), then it can be shown ^ that 
once again n^^av = {ns,A + ?^s,s)/2. 
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If the quasiparticle conductivity is frequency- 
dependent, then the analogous scalar results of 
the previous section continue to hold. For exam- 
ple, if aqp^A{>jj) = tTgp,_B(w), then the extra spec- 
tral weight due to the superfluid inhomogeneity is 
again given by eq. (5) in the weak-inhomogeneity 
regime. 

In the case of the stripe geometry, where the 
principal axes of the conductivity tensor take ei- 
ther of two perpendicular orientations with equal 
probability, n^^e is given by the duality result 

ns,e = ^/ns,Ans,B- (7) 

This form allows the extra spectral weight to be 
evaluated straightforwardly, given aqp{u!), with- 
out making the small fluctuation approximation. 

3. Numerical Example 

As a simple example, we consider a supercon- 
ducting layer in which the conductivity has one 
of two possible values, (Ta(w) or i7b(w) with equal 
probability, and we assume aA.si^^) = o'qp{uj) + 

2 

q "si^.-B where we take n, 4 > n, r. This model 
would be suitable either for a layer with static 
scalar disorder, or for a model of stripe domains, 
as discussed above. We also write aqp = cto/(1 — 
iojTqp), with CTo = Tiqpq^Tqplm* . We also assume 
Uqp — aT for T < Tc ot Uqp = aTc for T > Tc, 
Tqp^ = (3T, fis^A = 7"-s,o, "-s,B = 7^^ri5,o, where 
"-s,o(r) = ns,o(0)\/l - 2aT/nsfi{0), and 7 is a 
parameter describing the superfluid inhomogene- 
ity (this form ensures that Ug^av = ^Jns^Ans,B = 
Usfi). Our choice of temperature dependence for 
nsfi{T) ensures that Ugfi (i) decreases linearly 
with increasing temperature T at small T, as ob- 
served experimentally and (ii) vanishes at a 
critical temperature Tc as V^c — T. The ingre- 
dients of this model are very similar to those of 
Ref. and have a straightforward interpreta- 
tion. First, Us^qp should be proportional to T in 
a gapless d-wave superconductor and hence 
Us^e should be depleted by the same amount and 
fall off linearly in T . We also assume that Us^a 
and ria^B individually are linear in T . The form 
1/Tqp = j3T, where (3 is another constant has 
been observed for nodal quasiparticles in the su- 
perconducting states [p^. 



We compute Ue using Brug'geman effective- 
medium approximation (EMA) (5||9(] , which gives 

{a A - (Je)/{(yA + CTe) + [(7 B " CTe)/ {(T B + (Je) = 0. 

The solution to this equation is simply = 
^JoA^B- Fig. 1 shows the resulting ae.,\(y>^T) 
for several frequencies ranging from 0.2 to 0.8 
THz, the range measured in Ref. ||]. The pa- 
rameters Tc, (To J ^s,o(0), a and /3 were taken 
from Ref.^, and we assumed 7 = 8. Also 
shown are J^"""' ae,i{Lu,T)dLu for tj,„i„/(27r) — 

0. 2.Hz and d^^ax/i^ir) = O.STHz. Finally, we 
plot (TqpA{uJ,T) for these frequencies, as well as 

Cgp^dw. Clearly, (Te,i('^) is considerably 
increased beyond the quasiparticle contribution, 
because of spatial fluctuations in the superfluid 
density. 

4. Discussion 

We have shown that a superconducting layer 
with an inhomogeneous superfluid density will 
have an extra absorption not present in a homoge- 
neous superconductor. The frequency integral of 
(Te,i(w) associated with this absorption is propor- 
tional to the superfluid density, in general agree- 
ment with experiments [0J|]. We have also shown 
that this inhomogeneity, and hence the extra ab- 
sorption seen in experiments, can arise from a 
stripe domain structure, as well as from random 
but isotropic disorder. Thus, we speculate that 
no such extra absorption should be observed in 
a high-Tc superconductor unless one of these two 
types of inhomogeneities are present (beyond that 
expected purely from nodal quasiparticle absorp- 
tion). 

5. Acknov^rledgments. 

This work has been supported by NSF Grant 
DMR01-04987, and by the U.-S. /Israel Binational 
Science Foundation. 

REFERENCES 

1. J. Corson et al, Nature (London) 398, 221 
(1999). 

2. S.-F. Lee et al, Phys. Rev. Lett. 77, 735 
(1996). 



4 



a [|a.n m] 

e,l 



1.2 

1 

. 8 
. 6 
0.4 
0.2 



-1 



(a) 




20 40 60 80 100 T, K 



Z [[iD. m] THz 

e,l (b) 




20 40 60 80 100 



T,K 



a [\\D. m] 

qp.i 



(c) 



Z m] THz 

qp.i (d) 



0.7 




20 40 60 80 100 T, K 




20 40 60 80 100 



T,K 



Figure 1. (a) (j,,^i{uj.T), for oj/{2tt) = 0.2 THz (solid line), 0.4 THz (dotted line), and 0.8 THz 
(dashed line), for the model inhoniogcneous superconductor described in the text. Also plotted are (b) 
Se,i = Jl^^°^ (7e,iioj)duj, (c) CTgj,,i(w,T), and (d) S,p_i = J^^"-^ agp^i{uj)dLO, where a;mj„/(27r) = 0.2Thz 
and LOmaxJl^T^) = O.STHz. 
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